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Abstract 


We study scalar perturbations of nonlinear charged Lifshitz black branes with hyperscaling vi¬ 
olating factor, and we find numerically the quasinormal modes for scalar fields. Then, we study 
the stability of these black branes under massive and massless scalar field perturbations. Also, we 
consider different values of the dynamical exponent, the nonlinear exponent and the hyperscaling 
violating exponent. 
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I. INTRODUCTION 


Lifshitz spacetimes have received considerable attention from the condensed matter point 
of view due to the AdS/CFT correspondence, i.e., the searching for gravity duals of Lifshitz 
fixed points for condensed matter physics and quantum chromodynamics [l|. From the 
quantum field theory point of view, there are many invariant scale theories of interest when 
studying such critical points. Such theories exhibit the anisotropic scale invariance t —)■ 

X x^i with 2 : 7^ 1, where z is the relative scale dimension of time and space, and these are 
of particular interest in studies of critical exponent theory and phase transitions. Systems 
with such behavior appear, for instance, in the description of strongly correlated electrons. 
The importance of possessing a tool to study strongly correlated condensed matter systems 
is beyond question, and consequently much attention has focused on this area in recent 
years. In this sense, Lifshitz holographic superconductivity has been a topic of numerous 
studies and interesting properties are found when one generalizes the gauge/gravity duality 


to non-relativistic situations Hm. 

The Lifshitz spacetimes are described by the metrics 


„2z 


£2; 


ds^ = -^-Zde + 


£2 


( 1 ) 
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where x represents a. D — 2 dimensional spatial vector, D is the spacetime dimension and ^ 
denotes the length scale in this geometry. If 2; = 1, the spacetime is the usual anti-de Sitter 
metric in Poincare coordinates. Furthermore, all scalar curvature invariants are constant 
and these spacetimes have a null curvature singularity at r ^ 0 for 2; 7^ 1, which can be 
seen by computing the tidal forces between infalling particles. This singularity is reached 
in hnite proper time by infalling observers, so the spacetime is geodesically incomplete 12 |. 
The metrics of Lifshitz black holes asymptotically have the form ([1]); however, obtaining 
analytical solutions does not seem to be a trivial task, and therefore constructing finite tem¬ 
perature gravity duals requires the introduction of strange matter content with a theoretical 
motivation that is not clear. Another way of hnding such a Lifshitz black hole solution is by 
considering carefully-tuned higher-curvature modihcations to the Hilbert-Einstein action, as 
in new massive gravity (NMG) in 3-dimensions or corrections to g eneral relativity. Some 


Lifshitz black holes solutions have been found in references 


13l-ll9| . Thermodynamically, 


it is difficult to compute conserved quantities for Lifshitz black holes; however, progress 
was made on the computation of mass and related thermodynamic quantities by using the 


ADT method [3, l2l| and 
Euclidean action approach 


he off-shell extension of the ADT formalism 2^ as well as the 


m 


24i\ . Also, phase transitions between Lifshitz black holes and 


other configurations with different asymptotes have been studied in 2^. However, due to 
their different asymptotes these phases transitions do not occur. 

A generalization of the above metric is given by 


r.2z 


ds^ = r D I --—dr — -dr^ + -r^dx 




P 


( 2 ) 


which, besides having an anisotropic scaling as the Lifshitz metric, have an overall hyper¬ 
scaling violating factor with hyperscaling exponent 0, thus, this line element is conformally 
related to the Lifshitz metric. This space-time is important in the study of the dual field 


theories with hypersca^ 
have been found in 28 


i ng violation 

4Qi. 
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27| . Lifshitz black holes with hyperscaling violation 


In this work, we study scalar perturbations of nonlinear charged Lifshitz black branes 
with hyperscaling violation. The matter is parameterized by scalar fields minimally coupled 
to gravity. Then, we obtain numerically the quasinormal frequencies (QNFs) for scalar fields. 
We focus our study in the influence of the dynamical exponent, the nonlinear exponent and 
the hyperscaling exponent in the stability. 
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The study of the QNFs 3lM36[| gives information about the stability of black holes under 


matter helds that evolve perturbatively in their exterior region, without backreacting on the 
metric. In general, the oscillation frequencies are complex, where the real part represents 
the oscillation frequency and the imaginary part describes the rate at which this oscilla¬ 
tion is damped, with the stability of the black hole being guaranteed if the imaginary part 
is negative. The QNFs are independent of the initial conditions and depend only on the 
parameters of the black hole (mass, charge and angular momentum) and the fundamental 
constants (Newton constant and cosmological constant) that describe a black hole, just like 
the parameters that dehne the test held. On the other hand, the QNFs determine how 
fast a thermal state in the boundary theory will reach thermal equilibrium according to 
the AdS/CFT correspondence [^, where the relaxation time of a thermal state is propor¬ 
tional to the inverse of the imaginary part of the QNFs of the dual gravity background, 
which was established due to the QNFs of the black hole being related to the poles of the 
retarded correlation function of the corresponding perturbations of the dual conformal field 


theory 38|]. Fermions on a Lifshitz background were studied in 39| by using the fermionic 


Green’s function in 4-dimensional Lifshitz spacetime with z = 2; the authors considered 
a non-relativistic (mixed) boundary condition for fermions and showed that the spectrum 
has a flat band. Also, the Dirac quasinormal modes (QNMs) for a 4-dimensional Lifshitz 
black hole were studied in j^|. Generally, the Lifshitz black holes are stable under scalar 


pertur 


m 
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Dations, quasinormal modes under scalar held perturbations have been studied in 
- l48| and electromagnetic quasinormal modes in 49(]. Moreover, it is was stablished 


that for d> z + 1 , aX zero moment^ 


he modes are non-overdamped, whereas for d < z + 1 


the system is always overdamped The QNFs have been calculated by means of nu¬ 


merical and analytical techniques, some remarkably numerical methods are: the Mashhoon 
method, Ghandrasekhar-Detweiler, WKB method, Frobenius method, method of continued 
fractions, Nollert, asymptotic iteration method (AIM) and improved AIM among others. In 
the context of black hole thermodynamics, QNMs allow the quantum area spectrum of the 
black hole horizon to be studied as well as the mass and the entropy spectrum. 


The paper is organized as follows. In Sec. [TTlwe give a brief review of nonlinear charged 
Lifshitz black branes with hyperscaling violation that we will consider as background. In 
Sec. im we calculate the QNFs of scalar perturbations numerically by using the improved 
AIM. Finally, our conclusions are in Sec. IIVI 
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II. NONLINEAR CHARGED LIFSHITZ BLACK BRANES WITH HYPERSCAL¬ 
ING VIOLATION 


The nonlinear charged Lifshitz black brane that we consider is a solution of the Einstein- 
dilaton gravit y in the presence of a linear and a nonlinear electromagnetic held, this solution 


was found in 




29| . The action is given by 

167r Im + V^c^) - )‘^ , (3) 

where R is the Ricci scalar on the manifold AI, 0 is the dilatonic held, Ai and A 2 are free 
parameters of the model, F is the Maxwell invariant of the electromagnetic held 
where is the electromagnetic potential and = d^^By^ is a linear electromagnetic held, 
where is the electromagnetic potential. The following metric is solution of the theory 
dehned by the action ([3]), and it represents a black brane solution with hyperscaling violating 
factor 


ds^ = r 


2a 


r f{r)dt + 




D-2 


with 


f{r) = 1 - 


r‘^f{r 
M 


+ ^ dxl 


2=1 


^z+D-2-6 


+ 


.T+D-2+z-e ’ 


(4) 


(5) 


where a = been used, and 9 is the hyperscaling exponent, also 

n 2 s ^ _J_o_ (2a'^Y 

^ 4{D-2-e)r ^ ’ 

where M is an integration constant related to the mass of the black brane and q 2 is an 
integration constant related to its electric charge. The solutions are not valid for a = —1. 
To have /(r) ^ 1 when r ^ 00 , the following inequalities must be satished 


z + D -2-9 > 0 , T + D-2 + z-9>0, 


(7) 


however, it can be shown that T > 0 29(]. The gauge and dilatonic helds are given by 


Fr, = 


Hrt = qiro^ r 


{D-2+Z-1) 


^2{D-2-e){-T^+z-l) 


(j){r) = In 


5 


( 8 ) 












thus, for a real dilatonic field we must have {D — 2 — 9) {—^Z 2 +2^ — 1) >0. Moreover, 
since z + D — 2 — 9>Q, z has to be larger than 1 too. It is worth to mention that the 
condition for having a black hole is 

/ y+^+D-2-0^ . pg2s 

\T + z + D-2-e) - \z + D-2-e 

The temperature and entropy of the solution are given by 

r = ^ ((^ + £> - 2 - e)Ti - . s = . (lo) 


z+D-2-e 


(9) 


where Vh denotes the event horizon, 
detail in 29 1 . 


The study of the thermodynamics was performed in 


III. QUASINORMAL MODES 


The QNMs of scalar perturbations in the background of a D-dimensional nonlinear 
charged Lifshitz black brane are given by the scalar field solution of the Klein-Gordon 
equation with suitable boundary conditions for a black brane geometry. This means there 
are only ingoing waves on the event horizon and we consider that the scalar field vanishes 
at spatial inhnity, known as Dirichlet boundary condition. The Klein-Gordon equation for 
a scalar held minimally coupled to curvature is 

, ( 11 ) 

where m is the mass of the scalar held <p. Now, by means of the following ansatz 

(p = , ( 12 ) 


where T is a spatial vector in H — 2 dimensions, and —k‘^ is the eigenvalue of the Laplacian 
in the hat base submanifold. The Klein-Gordon equation reduces to 


1 d 




U 


K 


J ^ ^2 q -|-22^ (^2 d -|-2 


— ) R{r) = 0 , 


(13) 


where (3 = {a + 1)D + z — ?> has been dehned. Now, dehning R{r) as 


( 14 ) 
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where n = and by using the tortoise coordinate x given by 

dr 

= r-+if^r) ’ 

the Klein-Gordon equation can be written as a one-dimensional Schrodinger equation 

d‘^F{x) 


( 15 ) 


dx"^ 


V{r)F{x) = -uj^F{x ), 


(16) 


with an effective potential V (r) given by 

V{r) = + {D — 2)(1 -f a) ^((-D — 2)(1 -t- a) -f 2z) /(r) -|- 2r^ ) r 

(17) 

that diverges at spatial infinity. It is worth to mention that is not trivial to find analytical 
solutions to Eq. m In the next section, we will perform numerical studies by using 
the improved AIM Ml], which is an improved version of the method proposed in references 


52 


53l | and it has been app 


geometries, see for instance 


lec 


40 


successful in the context of QNMs for different black holes 


51 


54-5d. 


1. Stability analysis 


Following the argument used in 57||, adapted to Lifshitz geometries with hyperscaling 
violation, we can verify when the imaginary part of the quasinormal frequency uj is always 
negative. By using outgoing Eddington-Filkenstein coordinates v = t + x, metric (jl]) can be 
transformed to 


D-2 


ds^ = I r‘^^f{r)dv'^ + 2r^ ^dvdr + ^ dxf 

Now, taking as ansatz 


i=l 




with n = 2 )(i+q!) ^ Klein-Gordon equation yields 


A. 

dr 


[F+^f(^r)ijj'{r)) — 2iu‘ip'{r) — V{r)'ip{r) = 0 , 


where 


(18) 


(19) 


( 20 ) 


V (r) = nr^f{r) -|- n(n -|- z)r^ ^ f{r) -|- ^ ^ 


( 21 ) 
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Notice that n > — | according to the inequalities ([7]). Then, multiplying equation (|20|) by -0* 
and performing integrations by parts, and using Dirichlet boundary condition for the scalar 
held at spatial inhnity, one can obtain 




‘+7(r) 


dijj 

dr 


+ V{r) 


OJ? |V^(r = rfe)|^ 
Imioj) 


( 22 ) 


thus, the stability of the black brane under scalar held perturbations is guaranteed for a 
strictly positive potential V(r) outside the horizon, because in this case, equation fl22|) is 
satished only for Imioj) < 0. Notice that the potential fl2T|) is positive for n > 0 or a > —1, 
what guaranties the stability of the black brane solution. In this work we focus our attention 
to this case (a > —1). 


2. Numerical analysis 


In order to implement the improved AIM we make the following change of variable y = 
1 — r/j/r to equation flT3|) . where r^ denotes the location of the event horizon, thus, in this 
coordinate, the event horizon is located ai y = 0 and the spatial inhnity at y = 1. Then, 
the Klein-Gordon equation becomes 

d^R ( 2a-B f{y)\ m 

dy'^ V 1-1/ fiy)J dy 


a;^(l — yY^ ^ 


2^2a 




i? = 0, 


rl^f[yY 

in this equation /(?/) refers to the function /(r) evaluated at r = that is 


(23) 


f{y) 


M (1 - y)^+D-‘2-e _ yy+D-2+z-e 

z-\-D —2 —6 r+Z)— 2-\-z—6 

' h ' h 


(24) 


and f'iy) = Now, we must consider the behavior of the scalar held on the event 

horizon and at spatial inhnity. First, notice that the tortoise coordinate is given in terms of 
the y coordinate by 


dr 

dx = — , 
r^+^/(r) 


{l-yY ^dy 

rlfiy) 


(25) 


• Event horizon 


In the limit ?/ —)■ 0 (r ^ r/j) the function f{y) tends to f{y) = f\y)y + C>(|/^), where 
0{y’^) denotes terms of order y'^ and higher which can be neglected when y —)■ 0. 
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(26) 


Therefore, the tortoise coordinate is given explicitly by 

ln{y) 


X = 


/'(0)r 


where /'(O) = and the effective potential V{x) tends to zero in this limit; 

thus, equation (fT6|) reduces to 


(PF{x) 

dx"^ 


—uj^F{x ), 


(27) 


and its solution is 

F{x) = Cie-^^ + . (28) 


Imposing as boundary condition that only ingoing waves exist on the event horizon, 
we must set (72 = 0. Therefore, the solution near the horizon is given by 

_ itj 

F(x) = ^ y . (29) 


Spatial inhnity 

On the other hand, when ?/ —)• 1 (r —)■ cxo), the tortoise coordinate is given by 

( 1 - 2 /)" 


X = 


zrf 


(30) 


In the following we will consider two cases: the hrst case corresponds to — 1 < a < 0, 
and the second case corresponds to a > 0 and m = 0. In both cases the effective 
potential tends to: 


V{x) = 


6 


z2x‘i 


(31) 


where 5 = n{n + z) = 1(1 + a){D — 2){2z + (H — 2)(1 + a)). Then, equation flT^ for 


a; —)■ 0 becomes 


whose solution is 


d^F{ 


X) 


dx"^ z^x"^ 


F{x) = 0, 


F{x) = ^ 


(32) 

(33) 


Notice that the effective potential asymptotically tends to +cx) or to —cx) depending 
if 6 is positive or negative, respectively. In this work we focus on 5 > 0, and imposing 
Dirichlet boundary condition, that is, to have a null held at spatial inhnity, we must 
set Hi = 0. Therefore the solution becomes 


9 


rsj 


(34) 














Notice that, due to ([7]), n-\- z > zj^. So, requiring b = n{n + z) > 0 implies n > 0, or 
a > — 1 . 


Thus, taking into account these behaviors we dehne 


R (y) = y-^(l - + V + , 


(35) 


as ansatz. Then, by inserting these helds in Eq. (|2^ we obtain the homogeneous linear 
second-order differential equation for the function x{^) 


where 

My) = 


X" = My)x' + So{y)x , 


-Mhf'Wyi^ -h 2a - D{1 a) -1- z\/l + 46/z'^) 


(36) 


r?f'{0)y{l — y)f{y) 

+2i{l - y)uj)f{y) - rlf{0)y{l - y)f{y )), 

_^ . .. f2r;/'(0)W(l - »)"<•+“> - rl(l - y)^ 


(37) 


2rf+'>/'(0)%'^(l - !/)2 Ii+”)/(!/)2 
(rf/'(0)V(-2^ + {D- 2)z{l + a)(l + + 45 /^ 2 )) - 2ir^/'(0)(l - y) 

2)(1 + a)y + z^/l + AS/zM^ + 2(1 - 1 /)V)/(?/)^ - rlf{0)yf{y) 

(2r-ft/'(0)|/(r2“+W + K^{1 - yMM + r^(l - 2/)^+2“(r^/'(0)z(l + ^l + Ab/zM 

+2i{l-y)u)f\y))). (38) 

Then, in order to implement the improved AIM it is necessary to differentiate Eq. 
n times with respect to ?/, which yields the following equation: 


x""^ = My)x' + Sn{y)x , 


(39) 


where 

My) = Mi{y) + SnMy) + My)Mi{y ), (40) 

Sn{y) = s'n-Ay) + My)Mi{y) ■ (4i) 

Then, by expanding the and in a Taylor series around some point y, at which the 
improved AIM is performed yields 

OO 

Mv) = J2^n{y-vT , (42) 

i=0 
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( 43 ) 


i=0 

where the and are the Taylor coefficients of Xn{v) ^-nd Snirf), respectively, and by 
replacing the above expansions in Eqs. fHU]) and flTT]) the following set of recursion relations 
for the coefficients is obtained: 

i 

4 = (»+ i)4y + 4-1 + E =54-'i. («) 

A:=0 


if„ = (i + i)<;;.«+E‘'' 


k i—k 
O^n-1 


(45) 


fc =0 


In this manne 

rh 

t. 

le 

the AIM in 

51 

5 

541 


54| . and the quantization condition, which is equivalent to imposing a 


termination to the number of iterations, is given by 

ificy, - = 0. 


(46) 


We solve this equation numerically and we choose different values for the parameters. Thus, 
without loss of generality, we choose the following values D = 4, M = 4, z = 3, ro = 1, 
m = 0 and k = 0 in Tables HI UTI and HTTl In Table HI we show some lowest QNFs, for massless 
scalar field with s = 1.5, q 2 = ^ and different values of 6. Then, in Table m we show some 
lowest QNFs, for 6 = 1, q 2 = 1 and different values of s and in Table HTTI we show QNFs with 
9 = 1, s = 2 and different values of q 2 - In Table HVl we show QNFs for H = 4, M = 4, 6* = 1, 
s = 2, tq = 1, q 2 = I, m = 0, K = 0 and different values of 2 ;. We observe that in all cases 
analyzed the QNFs have an imaginary part that is negative, which ensures the stability of 
nonlinear charged Lifshitz black branes with hyperscaling violation under massless scalar 
perturbations. Then, in Table IVl we show some lowest QNFs for massive scalar helds, for 
s = 2, 6 = 1, K = 1, and different values of m. Finally, in Table |Vl] we show fundamental 
QNFs for massless scalar held for D = A and D = 5, with different values of the dynamical 
and hyperscaling violating exponents. Note that the Klein-Gordon equation depends on the 
combination D — 6, this is the reason why, for instance, the same QNFs are obtained for 
D = 9 = 0.5 and D = 5, 9 = 1.5 in Table IVTl 
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Table I. Some quasinormal frequencies for D = 4, M = 4, s = 1.5, rg = 1, ^2 = 1; 2 : = 3, m = 0, 


K = 0 and different values of 0. 


n 9 = —1 


= -0.5 


0 = 0 


9 = 0.5 


9 = 1 


= 1.5 


0 -10.71750i -8.77694i -7.57904i -6.57065i -5.45750i -2.85448i 

1 -12.82510i -13.22240i -12.62600i -11.65420i -10.15690i -5.56952i 

2 -21.81980i -19.50890i -18.14820i -16.81610i -14.83420i -8.26014i 

3 -25.13070 -24.875901 -23.579701 -21.992201 -19.516701 -10.944601 


Table II. Quasinormal frequencies for H = 4, M = 4, z = 3, 0 = 1, rg = 1, 92 = Ij w, = 0, k = 0 
and different values of s. 


n 

s = 1.5 

s = 2 

s = 2.1 

0 

-5.45750i 

-3.16787i 

-2.22508i 

1 

-10.15690i 

-6.05288i 

-4.30905i 

2 

-14.83420i 

-8.87462i 

-6.34705i 

3 

-19.51670i 

-11.68160i 

-8.36704i 


IV. CONCLUDING COMMENTS 

In this work we have calculated numerically the QNFs of scalar held perturbations of 
nonlinear charged Lifshitz black branes with hyperscaling violation by imposing suitable 
boundary conditions on the event horizon and at spatial inhnity. The scalar held is con¬ 
sidered as a mere test held, without backreaction over the spacetime itself. In this study 
we have considered the case a > — 1, additionally, for a > 0 we have restricted to massless 
scalar held m = 0. Then, we have studied the stability of these black branes under massive 
and massless scalar held perturbations. In general, our results show that the QNFs have a 
negative imaginary part. Therefore, the black brane is stable under massive and massless 
scalar held perturbations. Also, we can see that there is a limit on the dynamical exponent 
z above which the system is always overdamped for a given dimension, and the hyperscal¬ 
ing violating exponent shifts this limit. For instance, as we can observe in Table IVIl for 
D = 5 and 0 = —1 we have found that the system is non-overdamped for z = 3.5, but it is 
overdamped for z = 4. However, for D = 5 and 0 = 0 we have found that the system is non- 
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Table III. Quasinormal frequencies for L) = 4, M = 4, 2 : = 3, s = 2, 0 = 1, tq = 1 m = 0, k = 0 
and different values of q 2 - 

n q 2 = 0.1 q 2 = 0.5 ^2 = 1 

0 -6.66991Z -6.49942i -3.16787i 

1 -12.25930i -11.96360i -6.05288i 

2 -17.91470Z -17.48060i -8.87462i 

3 -23.57020Z -22.99920i -11.68160i 

Table IV. Quasinormal frequencies for D = 4, M = 4, 0 = 1, s = 2, ro = 1, ^2 = 1, = 0, k = 0 

and different values of 2 . _ 

n z = 3 z = 5 z = 8 

0 -3.16787Z -8.80358i -15.17150i 

1 -6.05288i -16.87330i -29.48100z 

2 -8.87462Z -24.90450i -43.77070z 

3 -11.68160i -32.93630i -58.06180i 

overdamped for = 2.5 and it is overdamped for z = 3, and for D = 5 and 6 = 0.5 we have 
found that the system is overdamped for = 2.5. It is worth to mention that the shift also 
depends on the dimension. On the other hand, when we increase the hyperscaling violating 
exponent the relaxation time of the dual thermal states increases (due to the absolute value 
of the imaginary part of the QNFs decreases) and when we increase the nonlinear exponent 
the relaxation time increases too, when the system is overdamped. 
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